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A recently introduced discrete formalism allows to solve the problem of time in quantum gravity 
in a relational manner. Quantum mechanics formulated with a relational time is not exactly unitary 
and implies a fundamental mechanism for decoherence of quantum states. The mechanism is strong 
enough to render the black hole information puzzle unobservable. 
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I. INTRODUCTION 

The "problem of time" in quantum gravity arises 
largely due to the presence of constraints in the theory, 
in particular the Hamiltonian constraint (see [J for a re- 
view). If one could eliminate the constraints almost all 
the conceptual problems with the problem of time can be 
eliminated. In particular, one can implement the quan- 
tization proposed by Page and Wootters in which one 
promotes all variables to quantum operators and chooses 
one of these quantum variables to be the "clock" . One 
then computes conditional probabilities for other vari- 
ables to take certain values when the "clock" variable is 
at a certain "time" . 

We have recently introduced an approach to general 
relativity that eliminates the constraints. The main idea 
is to approximate the theory by a discrete theory, much 
like it is done in particle physics when one approximates 
a field theory on the lattice. The novelty of our approach 
is that the discrete theories we construct are constraint- 
free, yet they can approximate general relativity. As a 
consequence, one can complete the Page- Wootters quan- 
tization of the discrete theories and introduce a relational 
time 1^. 

An immediate consequence of having a "quantum 
clock" variable in quantum mechanics is that the evolu- 
tion is not unitary 5, 6]. Both the clock and the system 
under study evolve unitarily and under the usual rules 
of quantum mechanics in terms of a fiducial background 
time n (we use the letter n to emphasize that we are 
working in a discrete formulation, though this is not cen- 
tral to the points discussed in this paper). This time is 
an inaccessible variable, we could only measure it if we 
had a perfectly classical clock. What we can measure 
are the dynamical variables of the problem, in particular 
t, the variable that describes the clock. This variable is 
represented by a quantum operator and it will have an 
expectation value and a dispersion. Upon evolution, the 
dispersion will increase. One can show that if one pre- 
pares the clock initially in a state in which t is highly 



peaked around a given value of the fiducial time n, the 
quantities under study (let us call them O) will evolve 
according to an approximate Heisenberg equation, but 
there will be corrective terms that imply that pure states 
evolve into mixed states. 

This kind of discussion is relevant to the black hole 
information puzzle The puzzle arises because one 
could have a pure state that undergoes gravitational col- 
lapse to form a black hole. The black hole will eventually 
evaporate leaving behind outgoing thermal Hawking ra- 
diation. It will therefore be in a mixed state, so somehow 
the initial pure state evolved into a mixed state. This is 
a problem in ordinary quantum mechanics. In quantum 
mechanics with a relational time, since pure states de- 
cohere anyway, one should evaluate if the rate of funda- 
mental decoherence is slower or faster than the process of 
black hole evaporation. If the state would have become 
totally mixed anyway due to fundamental decoherence 
by the time the black hole evaporates, then the puzzle is 
unobservable. 

We have carried out some preliminary calculations to 
yield an estimate of the fundamental decoherence and 
compare it to the black hole evaporation rate The 
rate of decoherence is related to "how classical" the clock 
one uses is. Therefore we addressed the question of what 
is the optimal clock that one can construct. Following ar- 
guments of Salecker and Wigner and further discussions 
by Ng and Amelino-Camelia 0, we established that the 
optimal clock one can find is a black hole. The accuracy 
with which a black hole can measure time is given by the 
frequency of its (quasi)normal modes, which scales as the 
inverse of the black hole mass. One would like therefore 
to have a black hole of small mass as a clock in order to 
have more accuracy. But there is a limit to this, if the 
mass is too small, the black hole-clock will evaporate too 
quickly for one to observe the physics of interest. There- 
fore one has an optimal accuracy one can achieve given a 
lapse of time that one wishes to measure. This limit on 
the accuracy is clearly only theoretical, in practice there 
will be other environmental factors that will affect the 
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accuracy of the clock. Even if one isolates the system 
from the environment, there are quantum uncertainties 
(for instance, in the position of the clock) that need to 
be taken into account. Here we will only concentrate on 
the above mentioned effect, since it is one of the effects 
of fundamental nature and can be viewed as an ultimate 
limit on the accuracy of a clock. 

The organization of this paper is as follows. In section 
II and III we review results (0, Q) on the evolution of 
conditional probabilities and the estimate of an optimal 
clock, respectively. In section IV we present a novel way 
of calculating explicitly the evolution of the conditional 
probability in closed form for the black hole information 
puzzle scenario. This generalizes previous work in which 
calculations were carried out to first non-trivial order in 
an expansion |^ . 



II. EVOLUTION OF THE CONDITIONAL 
PROBABILITY 

To give a more quantitative version of things, let us 
compute the time evolution of the conditional probability 
of measuring an observable O as a function of a real clock 
t (for further detail see 0]). 

We start by considering the conditional probability de- 
fined we mentioned above for an observable O to take a 
certain value o when the clock variable takes a value i, 

E„Tr(Ff(n)p) 

and one could have omitted the last Pt{n) using the 
cyclicity of the trace since we are assuming that Pt and 
Po commute. (If the observables (O and t) have con- 
tinuous spectra, the projectors in the above expression 
should be understood as integrated over the interval, i.e. 
Po{n) — J^^Po>(n)do' and similar for Pf.) The sum over 
all possible fiducial times n is due to the fact that we do 
not know for which value of n the variable t takes the 
value we want, and as the clock disperses there will be 
several values of n that correspond to t. We now intro- 
duce the hypothesis that the clock and the rest of the 
system interact weakly and write explicitly the evolution 
of the projectors in the step parameter n to get, 

V{o e Ao\t e At)p = (2) 
E„Tr (ul{n)Po{0)U2{n)ul{n)Pt{0)Ui{n)pi®P2) 

^ E„Tr(Pi(n)pi)Tr(p2) 

_ E„ Tr (u^^{n)PMU2{n)p2) Tr (ul{n)Pt{0)Ui{n)p, 

^ E„Tr(P,(n)pi)Tr(p2) 

Here we have assumed the system breaks into two subsys- 
tems, the clock (system 1) and the variables under study 
(system 2) and that the density matrix for the total sys- 
tem is a direct product p = pi x p2 and pi^2 are evolved in 
the fiducial time n by unitary evolution operators C/i,2- 



From this expression, using the cyclic property of the 
trace, we can identify the expressions of the density ma- 
trix evolved in relational time. We start by defining the 
probability that the measurement t corresponds to the 
value n. 



Vn{t) 



Tr 



(Pt{mi{T^)piUl{n) 
E„Tr(P,(n)pi) 



(3) 



and notice that ^„ 'Pn(i) = 1- 

We now define the evolution of the density matrix, 

P2{t) = Y,U2{n)p2Ul{7i)Vn{t), (4) 

n 

and noting that 

Tr {P2{t)) - 5]P„(i)Tr (pa) = Tr(p2) (5) 

n 

one can equate the conditional probability Q with the 
usual expression for a probability in quantum mechanics. 



V{o\t)p ^ 



_ TTiPo{0)p{t)) 



Mm) 



(6) 



where the projector is evaluated at i = since in the 
Schrodinger representation the operators do not evolve. 
From here on we drop the subscript 2 on the density 
matrix, since it is understood that we are discussing the 
density matrix of the system under study. 

It should be noted that all the sums in n, due to the 
assumption that the time variable is semiclassical, are 
only nontrivial in the small interval Ajn since outside of 
it, probabilities vanish. Something else to notice is that 
when we introduced the projectors, there was an integral 
over an interval. Therefore in the above expression for 
the evolution of the density matrix, this has to be taken 
into account. Since the interval At is arbitrary, one can 
consider the limit in which its width tends to zero, apply 
the mean value theorem in the integrals, and the interval 
in the numerator and denominator cancel out, yielding 
an expression for p{t) that is independent of the interval, 
and involves the non-integrated projector Pt(0). 

We have therefore ended with the standard probabil- 
ity expression with an "effective" density matrix in the 
Schrodinger picture given by p{t). In its definition, it is 
evident that unitarity is lost, since one ends up with a 
statistical mixture of states associated with different n's. 
We also notice that probabilities are conserved, as can be 
seen by taking © and integrating over x. We recall that 
p is not the normalized density matrix; the latter can be 
easily recovered dividing by the trace. 

We will assume that Vn{t) = f{t ~ t,nax{n)) with / a 
function that decays quite rapidly for values of t distant 
of the maximum ^max which depends on n. 

To manipulate expression Q more clearly, we will as- 
sume we are considering a finite region of evolution and 
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we are in the limit in which the number of steps in that 
region is very large. Wc denote the interval in the step 
variable n as going from zero to N with N a very large 
number. We define a new variable v = en with dimen- 
sions of time such that Ne = V with V a chosen finite 
value. We can then approximate expression by a con- 
tinuous expression, 

/ dvf{t-t„,aAv))p{v). (7) 

In this expression imax(t^) = tmax{n = v/e) and 

p{v)^U2 (n = vle)pUl {n ^v/e). (8) 

In all the above expressions, when we equate n — v/e it 
should be understood as n = Int(u/e), which coincide in 
the continuum limit. (Notice that strictly speaking we 
should write p{v/e) to keep the same functional form as 
for p{n), but we will drop the e to simplify the notation.) 
To simplify things further we will assume that we chose 
a physical variable as our clock that has a linear relation 
with V, i.e. tmax('f^) ~ V- In practice this is really not 
possible, there will be departures from this linearity and 
this is another effect that should be taken into account 
and will probably lead to further decoherence. 

III. AN OPTIMALLY CLASSICAL CLOCK 

We now need to make some assumptions about the 
clock. As we argued above, we use the ideas of Salecker- 
Wigner and Ng-Van Dam and Amelino-Camelia to 
suggest that the "most classical" clock one can build is 
a black hole. Briefly described, the argument goes as 
follows: if one considers an ideal isolated clock, it will lose 
accuracy as its wavefunction spreads. To try to diminish 
this spread, one can increase the mass of the clock. This 
process cannot continue indefinitely because eventually 
one will have enough mass to produce a black hole (trying 
to keep the density low by making a larger clock does 
not work since then one has to take into account that 
matter is elastic, etc). The black hole is therefore the 
most accurate clock from this point of view and it is also 
attractive as a fundamental clock given its fundamental 
nature (it is made of spacetime itself). The accuracy of 
such a clock is given by the quasinormal frequencies of 
the black hole, which scale as the inverse of the mass. 
That is, to have a more accurate black hole clock, one 
needs it to have a small mass. But a black hole of small 
mass evaporates due to Hawking radiation quickly. This 
creates a tension between these two requirements that 
leads to a formula that determines the best accuracy one 
can achieve in the measurement of a time Tmax, 



<5t ^tp^r,„axAp (9) 

where is Planck's time and from now on we choose 
units where h = c = 1. 



The reader might question why the spread of the wave- 
function of the clock limits its accuracy. After all, pre- 
sumably the clock is interacting with an environment 
which prevents the wavefunction from spreading. We ig- 
nore this effect, since this interaction is further source of 
inaccuracy in the clock (effects like this have been stud- 
ied in 10]) and wish to concentrate on the spread, which 
is an effect of fundamental nature, unrelated to the envi- 
ronment. 



IV. APPLICATION TO THE BLACK HOLE 
INFORMATION PUZZLE 

We need to make a quantum model of the black hole 
in order to study its decoherence. Here we will make a 
very primitive model. We assume the black hole hori- 
zon's area (or equivalently its energy) is quantized. This 
is usually assumed in quantum black hole studies and in 
particular it is predicted by loop quantum gravity. We 
choose a basis of states for the black hole labeled by the 
energy (area). The problem has some resemblance to the 
problem of an atom that is in an excited state and emits 
radiation to reach its fundamental state. If one considers 
the physical system under study to be the atom plus the 
radiation field, its evolution is unitary. One would expect 
a similar situation to hold for the black hole interacting 
with the gravitational and matter fields surrounding it. 
Here is where the paradox lies, since the evaporation pro- 
cess leads to loss of unitarity for the total system. Our 
model will include information about the black hole and 
the surrounding fields such that it starts its evolution in 
a pure state, and we will study its evolution according to 
the formalism developed in section II. We consider the 
system as described by a density matrix, 

P = J2pab\Mt)><Mt)\, (10) 

ab 

where 

\Mt)>=\E{t)+ea,Eo-E{t)> (11) 

and where the first entry in the bra (ket) represents the 
energy of the black hole at instant t, which changes with 
time in an adiabatic fashion, the constant Eq represents 
the mean value of the total energy of the system (which 
is conserved) and Eq — E{t) is the energy of the field at 
instant t. We consider the state to be a superposition of 
states of the black hole that differ in energy from E{t) 
by ta- To simplify the analysis we consider only a pair 
of levels of energy that are separated by an energy pro- 
portional to the temperature, as one would expect for 
an evaporating hole. Concretely, the characteristic fre- 
quency for this energy is given by 
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with Tmax the hfetime of the black hole (how long it takes 
to evaporate) and the subscript 12 denotes that it is the 
transition frequency between the two states of the sys- 
tem. Although this model sounds simple-minded it just 
underlies the robustness of the calculation: it just needs 
that the black hole have discrete energy levels character- 
ized by a separation determined by the temperature of 
the black hole. It is general enough to be implemented 
either assuming the Bekenstein spectrum of area or the 
spectrum stemming from loop quantum gravity pl| . We 
assume that we start with the black hole in a pure state 
which is a superposition of different energy eigenstates 
(there is no reason to assume that the black hole is ex- 
actly in an energy eigenstate, which would imply a sta- 
tionary state with no radiation being emitted; as soon 
as one takes into account the broadening of lines due to 
interaction one has to consider a superposition of states 
within the same broadened level with a time dependent 
separation with a similar behavior). Therefore the den- 
sity matrix has off-diagonal elements. 

We now compute the evolution of the two level model 
for the black hole using the formulas we developed in 
section II. We consider the off-diagonal matrix element 
of the density matrix in an energy eigen-basis, pi2. Its 
time evolution is given by, 



Pi2{t) = Pi2i0) J dvr^{t)e^pi^t J 0Ji2{T)dT 



(13) 



We can now compute the integral in the exponent, 
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<y5l2 = / UJi2{T)dT = -2 



(14) 



tp (Tmax - vY 



1/3 



[tpT; 



2 ]l/3 
max! 



To compute the evolution we need to provide a model 
for Vv{t)- We will assume it takes the form, 



1 



(15) 



where the function O is one if |w — i| < t/2 and zero 
otherwise, that is, a step function of width r centered at 
t. As we shall see the determination of the decoherence 
of the state does not depend on the particular form of 
the width as a function of one only needs to recall that 
the final width is given by the limit for the accuracy of 
the clock computed in equation (|5Jl. 



T{t) 



(16) 



t-T(i)/2 



To compute the integral, and evaluate it for the value at 
evaporation time Tmax we make the variable transforma- 



tion (Tniax — v)/tp = u^, and write. 



r(r„iax) V2(87r) V tp 



,2, / ^^u' 
duu t p exp ^ 

u \ 2 (Snf 



2/3N 



Pl2(T,nax) = -77^ 7 exp | { — | | X 



(17) 



with limits of integration U = ^ —T{T„^ax) / {"^tp). For 
a Solar sized black hole Tj^g^^/tp — (Mgun/Mp)'^, and 
therefore the integration limits are large. The integral 
can be evaluated in closed form in terms of Fresnel inte- 
grals, but it is more instructive to write the asymptotic 
form. The modulus of the integral behaves asymptoti- 
cally as \/ Msun/Mp. One therefore has an estimate for 
the modulus of the density matrix element behavior. 



|Pl2(rmax)| ^ |Pl2(0) 



M 



Ip 



Msu 



10-^«|pi2(0)|. (18) 



So for astrophysical black holes the puzzle is unobserv- 
able. One could still ask what is the situation for black 
holes that are smaller. We should recall that we have 
neglected several effects that further imply decoherence, 
so it is likely that the effect is larger than the estimate 
we present here. 



V. DISCUSSION 

The calculation we have carried out here differs from 
those of our previous papers. In made a first 

estimate of the decoherence in the context of the black 
hole information puzzle. In that first estimate we did 
not use an optimal clock and used a cruder model of the 
spectrum of the black hole (temperature independent). 
This calculation yielded that the quantum state did not 
decohere entirely by the time of evaporation, though it 
decohered in a significant amount. In we used an 
optimal clock and an improved model of the spectrum of 
the hole, but used only the first order expansion of the 
evolution equation for the state as valid throughout the 
whole evolution, namely. 



dpjt) 
dt 



(19) 



where a(t) encodes the information about the quantum 
fluctuations of the clock and it is related to the width T{t) 
by a{t) = ^r^/24. In we did an exphcit assumption 
for the form of r(t), defining T^(t) as the difference for 
the spread of the clock in the interval [0, Tmax] minus the 
spread in the interval [t, Tmax], that is 



tp 



2/3 



Tniax i 



2/3' 



(20) 
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One can then determine (T(i) in the expansion to be, 
cr{t)^tp/36(-^^Y\ (21) 

\ max / 

The calculation in Q yielded the remarkable result 
that it erased completely the information by the time 
evaporation occurs. In this paper we have integrated the 
full evolution equation and again one finds that there is 
a large level of decoherence by the time of evaporation. 
In the case of laboratory-like experiences our previous 
approach is justified as roughly an expansion in wi2T{t), 
where wi2 is taken as a the natural energy gap of the 
system. In order to see this we can integrate now IjlSf) 
for a time independent spectrum obtaining, 

p,.(0^2^e-^-*""^""'^^ (22) 

T(t) UJl2 

Expanding now in uji2T{t), 

Mt)\ = 2\pM0) r^^''f^ (23) 

-\Pi2mil-^u;!2r'{t)...) 

Noticing that if one is using the optimal clock for a 
total interval T, according to @ r(T) = iJ'^T'^l'^ and 
therefore, 

u^l2T^i^)=ty^u?,^T-'l\ (24) 

we can immediately compare with a similar calculation 
with a two level system, up to first order, obtained by 
using |19|) in |^ , where the level of fundamental decoher- 
ence is(l4|. 

Both expressions are in agreement and the formalism 
is consistent. The effect is too small to be observed 
in the lab, unless one can construct a system with a 
significant energy difference between the two levels. 
The most promising candidate systems would be given 
by systems of "Schrodinger cat" type. Bose-Einstein 
condensates could in some future provide a s yste m where 
the effect could be close to observability j^, [l3 ■ On the 
other hand, one could design an experiment where the 
effect is intentionally large, i.e. choosing a clock that 
does not behave semi-classically, as a proof of principle. 
In the case of a Black Hole, where the spectrum is 



[1] K. Kuchaf, "Time and interpretations of quantum 
gravity", in "Proceedings of the 4th Canadian confer- 



explicitly time dependent, the calculation in [g| can be 
also seen as an approximation of the exact result in the 
case of astrophysical black holes. Even though strictly 
speaking it is not a valid expansion for t — Tmax, one can 
show that it is a good approximation to the exact result 
provided r(t) << Tmax — t, region in which the phase in 
(|15|l can be expanded in series. This condition translate 

into imax << Tmax ~ tp^'^Tinax and therefore covers a 
good portion of the [OiTmax] interval for Mp/M << 1. 

Several caveats are in order. To begin with, it is clear 
that we have taken a very crude model for the black hole 
and a more detailed calculation is needed before one can 
completely write off the black hole information puzzle as 
an observable effect, but the present calculation provides 
good hope that the problem can indeed be solved. A 
realistic calculation seems somewhat beyond the state 
of the art. For instance, it is clear that the calculation 
should model quantum mechanically the black hole but 
also the fields it interacts with in a detailed way in the 
context of a theory of quantum gravity. There is also 
the issue that in these calculations we are neglecting the 
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are treating it as a "star that disappears" in the sense 
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in order to have a definite "evaporation time" to use in 
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will imply a change in the causal structure of space-time 
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loss should eventually be properly framed in a space-time 
context. This paper can only be viewed as a further step 
towards understanding how the imperfection in clocks 
can yield to loss of coherence in quantum mechanics. 

Summarizing, we have shown that unitarity in quan- 
tum mechanics only holds when describing the theory in 
terms of a perfect idealized clocks. If one uses realistic 
clocks loss of unitarity is introduced. We have estimated 
a minimum level of loss of unitarity based on constructing 
the most accurate clocks possible. The loss of unitarity 
is universal, affecting all physical phenomena. We have 
shown that although the effect is very small, it may be 
important enough to avoid the black hole information 
puzzle. 
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